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Abstract 


This paper studies the linear optimization problem subject to a system of 
bipolar fuzzy relation equations with the max-product composition oper- 
ator. Its feasible domain is briefly characterized by its lower and upper 
bound, and its consistency is considered. Also, some sufficient conditions 
are proposed to reduce the size of the search domain of the optimal solution 
to the problem. Under these conditions, some equations can be deleted to 
compute the minimum objective value. Some sufficient conditions are then 
proposed which under them, one of the optimal solutions of the problem is 
explicitly determined and the uniqueness conditions of the optimal solution 
are expressed. Moreover, a modified branch-and-bound method based on a 
value matrix is proposed to solve the reduced problem. A new algorithm is 
finally designed to solve the problem based on the conditions and modified 
branch-and-bound method. The algorithm is compared to the methods in 
other papers to show its efficiency. 


AMS subject classifications (2020): 90C70, 90C46, 90C26. 


Keywords: Bipolar Fuzzy Relation Equation; Linear Optimization; Max- 
Product Composition; Modified Branch-and-Bound Method. 


1 Introduction 


Sanchez [31] has firstly studied Fuzzy Relation Equations (FREs) and their 
associated problems. Then, many researchers investigated them from a the- 
oretical standpoint and in view of applications [26, 28, 32]. 
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Various approaches were designed to solve a system of FREs such as 
the algebraic method [21], the matrix pattern method [25], the universal 
algorithm [27], and the improved Lichun and Boxing’s method [39]. A com- 
prehensive review of their resolution methods has been expressed in [7] and 
references therein. An extended kind of FREs is a system of fuzzy Relation 
Inequalities (FRIs), which its solution set can be completely determined by 
finding its minimal solutions and maximum solution similar to FREs [16]. 
Their applications can be seen in the supply chain [38] and Peer-to-Peer data 
transmission network system [22, 36, 4]. 

The above FREs and FRIs are increasing in each of the variables. In some 
applications, for example, in an application of product public awareness in 
revenue management, we need variables with a bipolar characterization [9]. 
Some researchers introduced such bipolarity effects with the max-min compo- 
sition operator in the application [9]. Li and Jin [19] showed that checking the 
consistency of the system of bipolar max-min FREs is NP-complete. There- 
fore, the resolution of the linear optimization problem with constraints of 
bipolar FREs will be NP-hard. Yang [37] proposed a bipolar path approach 
to finding the complete solution set of the system. The characterization of the 
solvability of bipolar max-product FREs was investigated with the standard 
negation [5] and the product negation [6]. 


The optimization of objective functions with FRE constraints is an inter- 
esting research topic [8, 14, 17, 24, 29, 30, 33, 34, 35]. Fang and Li [8] firstly 
studied the linear programming problem provided to a system of the max- 
min FREs. They proposed an algorithm based on the branch-and-bound 
method with the jump-tracking technique for its resolution. The algorithm 
was extended to the problem with the max-product composition operator 
in [24]. Wu, Guu, and Liu, [35] improved their approach by designing an 
efficient procedure. In the procedure, the branch-and-bound method was ap- 
plied based on an upper bound for its optimal objective value. Hence, the 
procedure checks much fewer nodes to find the optimal solution with respect 
to Fang and Li’s approach. A necessary condition has been given to find an 
optimal solution to the problem with the max-min [34] and the max-product 
[14] composition operator. Then, the necessary condition was extended for 
fuzzy relation programming problem with the max-strict-t-norm composi- 
tion [33]. Three rules were presented to simplify the process of finding the 
optimal solution based on the necessary condition [34]. Li and Fang [17] in- 
vestigated the resolution and optimization of a system of FREs with Sup-T 
composition, and they generalized the most known results in literature and 
provided a unified framework for the resolution and optimization of the Sup- 
T equation. Recently, fuzzy relation programming was extended to optimize 
separable functions in [15]. Different kinds of fuzzy relation programming 
have appeared like the max-min fuzzy relation programming problem with 
the addition-min FRIs [4], linear optimization with the addition-min FRIs 
[12], and fuzzy relation lexicographic programming [40]. Recently, Zhou et 
al. [41] considered the problem of optimizing a nonlinear objective function 
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subject to a system of bipolar FREs with the max-Lukasiewicz triangular 
norm composition. They equivalently converted the problem to a 0-1 mixed 
nonlinear integer programming problem. It is NP-hard, and its resolution has 
high computational complexity. To increase the efficiency of algorithms, re- 
searchers focused on special classes of the bipolar fuzzy relation programming 
problems. Two important classes from the problem are linear [9, 20, 23, 3] 
and geometric [1, 2] with bipolar FRE constraints. For the first time, Freson, 
De Baets, and De Meyer [9] formulated the system of bipolar max-min FREs. 
They obtained the solution set of each of its equations. Using this point, they 
determined the solution set of a system of bipolar max-min FREs. This set 
can be characterized by a finite set of maximal and minimal solution pairs. 
They then studied the linear optimization problem subject to the system with 
a potential application of product public awareness in revenue management. 
They also found its optimal solutions based on the structure of the solution 
set of the system. Li and Liu [20] considered the problem with the max- 
Lukasiewicz t-norm. They converted the problem into a 0-1 integer linear 
optimization problem and solved it by using integer optimization techniques. 
However, the techniques may involve a high computational complexity. To 
overcome the point, Liu, Lur, and Wu [23] used the useful property that each 
component of an optimal solution can either be the corresponding component 
of lower or upper bound value and they proposed a simple value matrix with 
some simplification rules to reduce the dimensions of the original problem. 
To improve computational efficiency, some other rules were presented for 
more reduction of the dimensions of the problem with the max-parametric 
Hamacher composition operator in [3]. With regard to the above points, a 
modified branch-and-bound method was designed for the resolution of the 
problem in [3]. 


In this paper, the problem with the max-product operator is investigated. 
Its simplification procedures are completely different from the rules in [23, 3]. 
The motivations of this paper are the answers to the following questions: 
1- How do we can detect and remove the redundancy constraints in the bipo- 
lar FRE system? 
2- What are the sufficient conditions of optimality for a feasible solution to 
the original problem? 
3- What are the sufficient conditions for the uniqueness of the optimal solu- 
tion to the original problem? 
4- How can we design an efficient algorithm to solve the original problem 
with respect to the answers to questions 1 and 2? 
To find the answers to the above questions, two characteristic matrices are 
defined based on the components of lower and upper bound vector. Some 
sufficient conditions are presented to remove some equations. Moreover, some 
sufficient conditions are proposed which under them, one of the optimal so- 
lutions of the problem is determined. Then, the sufficient conditions for the 
uniqueness of the optimal solution are expressed. Furthermore, a modified 
branch-and-bound method based on a value matrix is designed to solve the 


410 Aliannezhadi and Abbasi Molai 


reduced problem. A new algorithm is proposed to solve the original problem 
based on the conditions and modified branch-and-bound method. 

The structure of this paper is organized as follows: Section 2 introduces 
the linear optimization problem subject to bipolar max-product FREs. We 
also investigate the characterizations of its feasible domain. In Section 3, 
the optimality conditions are presented for a feasible solution to the problem 
and some theorems are given to simplify and reduce the problem. Section 4 
proposes an algorithm to solve the problem. Some numerical examples are 
presented to illustrate the algorithm in Section 5. A comparative study is 
done with other methods to show the efficiency of the algorithm in Section 
6. Finally, conclusions are given in Section 7. 


2 Linear programming problem with bipolar 
max-product FREs 


This section is divided into two subsections. In the first subsection, the 
linear programming problem subject to bipolar FREs is formulated. The 
characterizations of its feasible domain are finally illustrated in the second 
subsection. 


2.1 Formulation of the problem 


Let AT = (a;;) and A~ = (a;;) be two m x n fuzzy relation matrices with 
0 < aj,,a;; < 1 foreachie I= {1,2,...,m} and j € J = {1,2,...,n}. 
Also, assume that b = (b1,---,bm)? € [0,1]™ and that ¢ = (c1,...,¢n) isa 
vector of cost coefficients, where c; > 0 for each 7 € J. In this paper, the 


following programming problem is considered: 


min Z(a“)= Disies: (1) 


st. AtogVA onr=b, (2) 
where x = (21,...,2n)’ € [0,1]" is the vector of decision variables to 
be determined and —a denotes the negation of x, that is, -# = (1 — 


@1,...,1—2%,)?. The operator of “o” represents the max-T,, composi- 
tion, where T; denotes the product operator. Moreover, S(At,A~,b) = 
{x € [0,1]” | At oxV A~ o nx = 5}, which consists of a set of solution vec- 


tors x € [0,1]” such that 


be “— _ . 
se max {aj,.23,0;;.(1—2,)} =b; forallie I. (3) 
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Problem (1)—(2) with the real cost coefficients can be converted to a prob- 
lem with nonnegative cost coefficients in a similar method to Subsection 3.3 
in [9]. Hence, without loss of generality, we assume that c; > 0 for each 
ged. 


2.2 The structure of the feasible domain of problem 


(1)-(2) 


A system of bipolar max-T, FREs At oa V A~ 072 = b is called consistent 
if its solution set, that is, S(At,A7,6), is nonempty. Otherwise, it is in- 
consistent. Now, we focus on the system of bipolar max-T, FREs (3), when 
S(At,A~,b) #0. 


Lemma 1. A vector x € [0,1]” is a solution for the system of bipolar max- 
T, FREs (3) if and only if max {a;;. 25, Qj; (1—2,;)} < }; for allie I 
and 7 € J, and for each i € J, there exists an index 73; € J such that 


max {a ii, D5, O55, jel —)} <0 


Proof. It is obvious. 


and b; ge 7€tland 7 € J, it is assumed that 
5 = 0, then we 


Remark 1. For any aj;,a i 


if aj; = 


define min{ by : 
Ges 
aj 


Lemma 2. For any a; and b; with i € I and Lg € J, the een of 


i: Ah 55 
_ 


max { aj; Lj, Q;;. aa 


min{ 2 a ; if vos £5, Qj; j;-(1 — 25) )} <b; = ag ies at least 


one of ‘two statements ca = 0ora,, = 0 holds. According A Remark 1, 
ee 


eae D5, Oi; ral 


oe, 1}. 


Proof. The proof will be divided into four cases as follows: 


Case 1. aj; 


# O and a;; £0. 
Case 2. a. =a,. =0. 


Case 3. aj, = 0 and a;; £0. 


Case 4. a”. #0 and ai; = 
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Case 1. It is eee that the inequality max {aj,.0,j,a;;.(1—aj)} < bj 
holds if and only if aj. jj Sb; and a;;.(1—2;) < bj. If b; A 0, then the 
recent inequalities can ee ee be peuiiiien as 1 — aS <azp< +e On 


ag as 
the other es according to the assumption, We feels 0 <a; < 1. The 
; < min{ 4 at 


Whenever 6; = 0, the inequalities imply x; < 0 ae x; > 1 which do not. hold 
for any x; € (0, 1]. 

Case 2. Since aj, = a;; = 0 and b; > 0, the inequality of max{a;j;. £5, A;3-(1— 
x;)} < b; holds if and only if x; € [0,1]. On the other hand, we have 


i ij 
quently, the result i also true in this CARs 


Case 3. Since ay = = 0, then min{ ah 
We now have the following subcases:. 1. b; = 0 and 2. bh # 0. In the 
first subcase, we have max{1 — ay = 1, that is, max{1 — = 


= 1. On the other hand, the inequal- 


ity max {aj Lj, Q;;. (1—2;)} < b; = 0 holds if and only if a;;(1 — 25) = 
0. This implies ca x; = 1. In the second subcase, the inequality of 
max { aj; £5,043. (1 — 25) i} < b; holds if and only if @;;.(1 ae < b;, which is 


a < #;. On the other Hand, we have 0 < a; <1. 


. b; 
xz; < min{ an? 


equivalent to max{1 — oe 


aj 


<1=min{4 at 


Remark 1. Hence, the eae? is dae in both subcases. 


Whenever $(At, A~,b) 4 0, the lower and upper bound on the solution 
set for the system of equations (3) can be determined using the following 
lemma. 


Lemma 3. [1] The vector of @ = (%1,...,%n)7 is the lower bound on the 
solution set of equations (3), 


ij 
3 € J. Also, the vector of #@ = (4#1,...,4n)7 is oe upper bound on 1 the 
solution set of equations (3), where £; = -min{ ue 


Furthermore, if there exists 7 € J such ine for pa te TI, a;; < 6;, then 
&; = 0 and if there exists j € J such that for each i € I, aj, < bj, then 
&; =1, that is, max 0 = 0 and min 0 = 1 are defined. 


Lemma 4. [1] Assume that $(A*,A~,b) #0 and that its lower and upper 
bound are % and @, respectively. If there exists j9 € J such that %;, = 24, 
then z;, = £;, = £;, for all e € S(At,A~,6). Also, the solution set of 
system (2) (or (3)) is the same to the following system: 


{emi 05 ,4;;-(1 — axj)} =b; foralliel—T, 


Lj 5 Ly < v5 for all j EJ— {jo}; Lig = Lig = Liqs 
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where @; and #, for all 7 € J are defined on the basis of system (2) (or (3)) 


and I = {i € I| max{aj, .a;,,a;;,-(1—2j,)} = bi} and T FO. 


We first consider the special case of b; = 0 for i € I. Since 
max 1 apties a;;.(1— a;)} is always nonnegative, the following inequality 
ety, a;;-(1 —a,;)} <b; can be converted into the equation 
+ 
ay 


max{a 
max {aj .a;, a;;.(1— a;)} = b;. We now express the following lemma. 

Lemma 5. Let 6; = 0 for i € J. A vector x is a solution for ith equation 
of the system (3) if and only if max{1 — %,0} < a; < min{ a 1} for each 


— > 3, 
aj; ai; 


ped. 


Proof. Considering Remark 1, the proof can be easily obtained by the proofs 
of Lemmas 1 and 2. 


Lemma 6. [1] Suppose that $(At, A~,b) #0, and that its lower and upper 
are &; and %;, respectively. Then the solution set of system (2) (or (3)) is 
the same to the following system: 


{eres a;;(1—2;)}=b; forallieI—Ip, 


Li; <2; for all j € J, 


where Ip = {i € I|b; = O}. Also, &; and 4, for all j € J, are defined on the 
basis of system (2) (or (3)). 


Furthermore, if 6; = 0 for each 7 € J, then we can easily obtain the 
solution set of equations (3) and an optimal solution for problem (1)-(2). 
These points are expressed in the following corollary. 


Corollary 1. If 6 = 0 in the constraint part of problem (1)-(2), then we 
have 
S(At,A~,b) = {x|%< a < &} and 2* = Z is an optimal solution to problem 


(1)-(2). 


Proof. It can be easily seen that S(At,A7~,b) = {x | & < x < &} with 
regard to Lemmas 3 and 5. The proof is completed by showing that «* = & 
is an optimal solution of problem (1)-(2). Since @ € S(At,A7~,b), c > 0, 
and the problem is minimization, then <% is an optimal solution of problem 
(1)-(2). O 


Without loss of generality, from now on, we will assume that <; < @;, for 
each j € J, and b; > 0, for each i € I. 


+ a 
Qj Oi 5 


=0. 


Remark 2. In this paper, if ay, = a;, = 0, then we define ae 
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I 
p, 


Re i. Xx, i 
0 ap a a a * 1/2 . 
Figure 1: Illustration of function f(x;) = maa{a;,. -@5,0;,.(1 — 23)} 


Lemma 7. For any a; and 6; with i € I and 7 € J, the equation of 


ij? a5 
Pag 
(1 as )} = b; has a solution if and only if wpa. <b < 
aj, +a; 


a,,}. Also, its solution set is determined with regard to the following 


max { aj, (oj Og. 
max{ay, 
cases: 


+ _ _ by r 
Case 1: If aj; < bj < aj, then S(aj;, say 05) = {in (2 1) }; 


ij? a5 


Case 2: If af, <b; <a, 


= = bi : 
jj, then S(a aj;,4;;, di) = { max (1 ee .0) t 


Case 3: If cane ub - <b; < min {aj, ay}, then 


Sasa) = { max (1 - = 0), min (4..1)}. 
Proof. For piven aed 


ij> 4; € [0,1], and &; > 0, the oS of the function of 
max {aj,.2;,a;;.(1— a3) sy and the solution set of S(aj.,a;,,bi) can be ob- 
served Lae Figure 1 and determined, easily. 


aig, Qi55 


Note that the vectors of & and & are only the lower and upper bound on the 
solution set of equations (3), respectively. They are not necessarily feasible 
solutions to the system of equations (3). Moreover, we have S(A*, A~,b) C 
{x | @< a < #}. Considering Lemmas 3 and 7, each equation in system (3) 
can be satisfied by #; or £;. In order to store these facts, the characteristic 
matrices of Q* and Q~ are defined below. 


Definition 1. Define two characteristic matrices Qt = (dij )mxn and Q7 = 
(Gi; )mxn such that for each i € I and j € J, we have 


ali fabty—bo oy goal if ay -%) =b, 
0 otherwise, Z 0 otherwise. 


Also, a series of index sets is defined as follows. 
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Definition 2. (é) [1] For the matrix QT, define 
T}(@) ={teT | xy = 2; and qi, = 1} and J (x) ={j€ J | aj = 2%; and qf = 1}. 
Also, for the matrix Q~, define 


IT, («) ={i el] | xj = %3 and q;, = 1} and J; (x) = {7 € J | xj = 3 and q;, = 1}, 
for each i € I and j € J. Furthermore, let [;(x) = I} (a) UT; (x), for each 
ged. 

- + _ 7tys) gt — J+ (5 — _ 7-/(%, - _ y-(% : 
(4) Let I = Ij (&), J = Jj" (@), let G7, = Tj (@), and let J; = J; (2), for 
each i€ J andj € J. 


Considering Lemmas 2 and 7 and the above concepts, we present a nec- 
essary and sufficient condition for the solution of system (3) (or (2)). 


Theorem 1. A vector x € [0,1]” is a solution for the system of bipolar max- 
T, FREs Atoa\V Aq ona = b if and only if <a#<4éand U L(x) =I. 
jes 


Proof. With regard to Lemma 1, we show that <a <%and UJ;(x) =I 
jes 

if and only if max {ajj-25,0;;- (1—2,;)} < 8, for all i € I and j € J, and for 

each 7 € I there exists an index j; € J such that 


max $0 2H azz. (1 — x;,)} = b;. Considering Lemmas 2 and 3, all the in- 


equalities max {Gy ty, Gael —2xj;)} < bj, for i € I and j € J, hold if and 


only if ¢ < a < &. We now focus on proving J I;(x) = I. The equality 
je J 
is true if and only if for each 7 € J, there exists an index 7; € J such that 


j 
j € J, we have YU J;(x) = J if and only if U (IF (2) UI; (x)) = I. Moreover, 
jed jes 
we have [J (UF (x) UJ; (x)) = J if and only if for each 7 € I there exists an 
jed 
index j; € J such that (aj, = %;, & Gi, = l)or (2;, = £5, & q;;, = 1), with re- 
gard to Definition 2. By Definition 1 and Lemma 7, it can easily be seen that 
for each i € I, there exists an index j; € J such that (x;, = @;, & qi, = 1) or 


max fags aj, 455, (1 —x;,)} = b;. Since I;(x) = I} (a) UI; (x), for each 


(xj, = %;, & qj), = 1) if and only if max | 48, isi, (l1- x,)} = b;. 


For each 7 € J, label the values of #; and x; with boolean variables of 
yj; and >y,;, respectively. The following theorem is used to determine the 
consistency of system AtoxV A7~ on7nr=b. 


Theorem 2. A system of bipolar max-T, FREs At ox V A7~ on7x = b 
is consistent if and only if its characteristic boolean formula C = (AC; is 
i€l 
well-defined and satisfiable, where C; = Vy; V V n7y;. 
jedy jedp 
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Proof. With regard to Definitions 1 and 2, the proof is similar to the proof 
of Theorem 2.5 in [18] 


Theorem 3. Let & and & be the lower and upper bound, respectively. Then 
LG cr ur, for alla € S(At,A™,b), forall je J. 


Exactly, we have J, (2)=IF (when x; = £;) or Ij(z)=I; (when x; = &;) or 
I;(z) =@ (when Z; < 2; < &;). 


Proof. The proof is obtained from Definitions 1 and 2. 


We are now ready to present some theorems to simplify problem (1)-—(2) 
in the next section. 


3 Some optimality sufficient conditions for problem 


(1-2) 


This section studies some optimality conditions for problem (1)—(2). One of 
its optimal solutions is found under the conditions. Moreover, some sufficient 
conditions are proposed to guarantee its uniqueness, and a closed form is 
proposed to determine it. In this section, it is assumed that S(A*,A7~,b) 4 0. 
The following lemma presents a useful property of the optimal solution of 
problem (1)-(2). 


Lemma 8. Consider the optimization problem of (1)—(2). Then there exists 
an optimal solution x* = (z7,...,2*)7 such that for each j € J either ri = 8; 


es. 
or v5 = Xj. 


Proof. The proof is similar to the proof of Lemma 4 in [20]. 


Some theorems are first presented to reduce the search domain of the 
optimal solution of problem (1)—(2). The dimensions of the matrices of Q* 
and Q~ can be reduced using these theorems. 


Theorem 4. Let T,, = {i € I\ {i:} | Jj? D Jf & Jp 2D Jz}, for i, € I. If 
for vector x, where @ <a < &, mas max{T,(a;;,%;),Tp(a;;,(1 — 2j))} = bi 
JE 


holds for i = 71, then for each i € T;,, we have 


mag maz{Ty(a4;,2)),Tp(a;;,(1— 25))} = by 


Proof. If for the vector x, where ¢ < x < &, the following equality 


mag maz{Tp(aj,,2;), Tp (ay, (1—<2;))} = b; holds for i = i, then there 


exist some j, € J such that max{T,(aj ,,,j,),Tp(a;,;,,(1—- vj,))} = ba. 
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Since these equalities hold only at the values of #;,, for j1 € Jj, (x), or %j,, 
for j1 € J; (x), then j; € J; (x)UJ;, (x) #0. Also, with regard to Definition 
2 J @) c Ji, and J;,(x) C Jj,. On the other hand, for each i € T;,, we 
have J;' D> Jj and J; D Jj, _ Therefore, for each i € T;,, J;*(x) C J;* and 
J; (@) C Jp pecauie je yout. (Q)Cae we Cd) and. Cd, 
for each i € T;,. Since j1 € Jj{(x) U Jj, (x), Jt (x) C Jj‘, and Jz (x) C Jj, 
for each i € T;,, then it is concluded that (j,; € Jj‘(#) and j; € Jj*) or 


(j1 € Jj (x) and ji, € J; ). So, we can write (xj, = #;, and qj = 1) 


or (xj, = &), and Gj, = 1), for each i € T;,. With regard to Definition 
Ly Tale af 2h) = = b; or Tp(a;;,,(1 — %5,)) = i, for each i € T;,. Thus, 
the equality max{T, (a Os Bix) Tp(a Gig) 1 — x;,))} = 6; holds true for each 
i€T;,. Since ¢ < x < & and max{T,(a By» i) Tp(aj;,,(1 — 23,))} = bi, for 
each i € T;j,, then the equality ve max{T,(az. a;;,0;),Tp(a;;,(1 — x5))} = bi 


holds true for each 2 € Tj,. 


The following corollary is a direct result of Theorem 4. 


Corollary 2. Under the conditions of Theorem 4, all the equations with 
numbers i € T;, can be removed from the matrices of Qt and Q~ once % and 
& have been obtained. 


In the next theorem, an equivalent system with system (2) is presented. 


Theorem 5. A binary vector u € {0,1}” induces a solution « = + Vu 
for the system of bipolar max-T, FREs At ox V A~ o 72 = b if and only if 
(Qt — Q-)u+ Q-en > em, where V=diag(# — %) and ex is a k-dimensional 
vector with the unit components. 


Proof. The proof is similar to the proof of Theorem 3 in [20]. 


Lemma 9. If there exists a pair i € I and k € J such that gf=q,, = 1, 


then the reduced system of (Qt’ - Q-')u at 6 el ee e;,,-1 and the system 
of (Q* — Q~)u+Q~en > em have the same solution set, where the matrices 


of Or ae and e’,,_, are obtained by removing the ith row of the matrices 
of QT, Q-, and e,, respectively. 


Proof. Let SQ, Q-,u) = {ue {0, | (Or i= Q™)ut Q- €n > @m} and 
S(Qh,Q™,u) = {ue {0,1} | (QY —Q™ ut Q en > ena}, where QT’, 

Q~, and e,,_1 are obtained from QT, a , and €;, by removing their row of 
i, respectively. We will now show S(Q*,Q7,u) = S(Q*’,Q~,u). Consider 
the ith equation of system (2). Since g,=qj, = 1 and uz € {0,1}, the ith in- 
equality of (Qt -—Q7 )u+Q~ en > em, or equivalently Qtu+Q7(e,—u) > em, 
is automatically satisfied. Therefore, S$(Qt+,Q-,u) = S(Q*’,Q~',u). 


Theorem 6. If there exists a pair 7 € J and k € J such that 7.1 
then we can remove the ith row in the computation of the minimum objective 
value. 
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Proof. It is obvious from Lemma 9 and Theorem 4 in [20]. 


Now, some sufficient conditions are presented to determine one of the 
optimal solutions of problem (1)-(2). First of all, we express the following 
definition. 


Definition 3. Define two index sets J; and I as follows: 
fy =| )ay and Ip=I\h. 
jet 
Corollary 3. If Iz = 0, then 7* = # is an optimal solution of problem 


(1)-(2). 


Proof. If UTI A then & is a feasible solution for the system of equations 
jes 
(2). Hence, we can assign #; to x; due to c; > 0, for each j € J. 


If the vector of & is a feasible solution of problem (1)-(2), then & is 
an optimal solution due to c; > 0, for each 7 € J. Otherwise, we have 
to construct a solution z* with elements xj = £; or x} = %j such that 
U L(a*) =I and the value of Z(x*) is the minimum objective value. With 
jed 


regard to these points, we present the following theorems. 


Theorem 7. If there exists an index k € {)} J;* such that the following 
i€ls 

conditions are satisfied 
1. forall je U a Ch(&_ — Ee) < cj(4; — Z;) and 

t€I2g 
ie i SUE 

jeJ 

ek 
then there exists an optimal solution z* = (2})j;<y for problem (1)-(2) as 
follows: 

& ifj=k, 

c= for all j € J. (4) 


Z; otherwise, 


Proof. We show that «* is an optimal solution of problem (1)—(2). It is 
enough to show that i) 2* € S(At,A7,b) and ii) Z(x*) < Z(x), for each 
ze S(At,A,b). 
i) With regard to the structure of vector x* and I; \Ij7 C U I, , the following 
je J 
zh 
equalities hold: 


Uae") =(Upug =(Unug. (5) 
JET jet JET 
jJ#k 
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On the other hand, k € () J;* implies that [jf D Ih =I\ U I; . Therefore, 
i€ Is jed 


ogee @ ieee ee (6) 
jed 
With regard to the expressions of (5) and (6), it is concluded that L I;(x*) = 
jed 
I. Since ¢ < «&* < ¢ and LU I;(a*) =I, the vector x* is a feasible solution 
jes 
for the system of equations (3) with regard to Theorem 1. 
ti) For each x € S(At,A~,b) and x 4 x*, we have xz = @, or there exists 
an index j € U J; \ {k} such that 2; = #;. 
i€l. 


we have 


If a, = &k, then there exists an index j’ € J\ {k} such that xj = @; due 

to « # x*. Hence, we have Z(a*) < Z(a*) + ej (4; — &j) < Z(ax). If there 

exists an index 7 € U J; \ {k} such that 2; = %;, then Z(x) > Z(x*) with 
1€I2 

regard to the condition 1 in Theorem 7. 


Note that the optimal solution introduced in the relation (4) is not nec- 
essarily unique with regard to condition 1 in Theorem 7 and c; > 0, for each 
j € J. Considering Theorem 7, some sufficient conditions are expressed in 
Lemma 10 that under them, problem (1)—(2) has a unique optimal solution 
and the optimal solution is explicitly determined. 


Lemma 10. If for each 7 € J\ U J;*, cj > 0 and there exists an index 


i€lo 
k€ () J;* such that the conditions 
1€ Ig 

1. forall j € U Jt \ {k}, cx(&e — Ex) < ¢;(#; — £;), and 
t€Ig 

2d, \i SU 
je 
i#k 


are satisfied, then the optimization problem of (1)—(2) has a unique optimal 
solution of «* = (x}) je, as relation (4). 


Proof. Since the assumptions of Theorem 7 hold, problem (1)—(2) has an 
optimal solution as the relation (4). We now show its uniqueness. By the 
assumptions of &; < #; and c; > 0, for each j € J, the condition 1 in Lemma 
10 implies that c; > 0, for each j € U J;*\{k}. This point implies that c; > 
t€l2 
0, for each j € J\ {k}, with regard to the assumption for all 7 € J\ U J;*, 
t€I2 
c; > 0. For each « € S(At,A~,b) and z # 2*, we have xz = &; or there 
exists an index j € LU J; \ {k} such that x; = @;. If zy = %,, then there 
i€I2 
exists an index 7’ € J \ {k} such that 2; = @j due to x # x*. Hence, we 
have Z(2*) < Z(x*) + ej (&j" — &)) < Z(x) due to the condition 1. If there 
exists j € U J;* \ {k} such that 7; = @;, then Z(x) > Z(x*) due to c; > 0, 
i€ To 
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for each j € J \ {k}. This shows the uniqueness of optimal solution «* for 
problem (1)-(2). 


The lemmas, theorems, and corollaries of this section are firstly applied to 
reduce the size of the original problem. If all of the components of the optimal 
solution of problem (1)—(2) were not determined, then we have to solve the 
reduced problem. To do this, we will explain the modified branch-and-bound 
method to find the rest of its components in the next section. 


4 A procedure for the resolution of problem (1)—(2) 


We first define a simple value matrix in this section. Applying the value ma- 
trix and some points, the branch-and-bound method with the jump-tracking 
technique is modified to solve problem (1)—(2). Finally, an algorithm is pro- 
posed for the resolution of problem (1)-(2). 


4.1 Modified branch-and-bound method 


In this subsection, we rewrite the objective function (1) as follows: 
n nm n n n 
Z(a) = aay = Soa; + Soi = PRA te i &;) + Sah 
j=l J=1 j=l j=l J=1 


Now, the optimal solutions of problem (1)—(2) are the same with the optimal 
solutions of the following problem: 


min Z(x) = decile; — £5), (7) 
st. AtogxVA on7nr=b, (8) 
a € (0, 1)”. (9) 


_ n 
It is necessary to recall that Z(a*) = Z(a*) + S>cj;é;. Therefore, we try 
j=l 


to find the optimal solution of problem (7)—(9). Since @ < x < #, for each 
xz € S(At,A~,b) and c > 0, we have Z(x) > 0 and Z(x) > YS cj#; for 
jes 
each x € S(At,A7~,b). As it was stated before, if ¢ € S(At, A~,b), then we 
can set % as an optimal solution. Otherwise, we have  ¢ S(At,A7~,b). In 
order to find the optimal solution of problem (1)—(2) (or (7)—(9)) according 
to Lemma 8, we have to set some @; instead of #;, 7 € J, in the vector & 
such that the new vector &* satisfy all equations (8)-(9) and minimize the 
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objective function (7). To do this, we first express a useful property of the 
objective function (7). 


Proposition 1. Let x = (x;)je7 be a vector in $(A*, A~,b) where x; = & 
and 2; = £, for two indices t and k such that t 4 k and vectors % and & 
are lower and upper bounds of system (8)—(9), respectively. Construct two 
new vectors 2’ = (x) jey and x” = (x7) je) such that x, = % and x = xj, 
for each j € J \ {t}, and af = &, and xj = xj, for each j € J \ {k}. If 


Cr (Lk = Ex) < Cr (ht = E+), then Zig") < Z(z"'). 


Proof. It is obvious. 


It is necessary to recall the importance of the form of problem (7)—(9) with 
respect to problem (1)—(2). With regard to Proposition 1, if c,@, > c&;, then 
we could not conclude that Z(2”) > Z(a’) or Z(x") < Z(a2’), but problem 
(7)—(9) gives us more information. If cx (% — Ze) < c¢(@: — Zz), then we have 
Z(x") < Z(x'). We use the useful property of the objective function (7) to 
present a modified branch-and-bound method to solve problem (1)—(2). First 
of all, we consider the following remark. 


Remark 3. Rearrange the rows of matrices Q* and Q~ such that the first 
|I>| rows in these matrices are the rows i € Ig. More precisely, transfer all 
rows i € Ig to the top |J5| of the rows of the matrices QT and Q-. 


According to Lemma 8, there exists an optimal solution x* = (27) jeJ 
such that either xj = %; or x; = 2%; for each 7 € J. Hence, it is concluded 
that Lait or Ij(2*)=I; , for each 7 € J, with regard to Theorem 3. 
Since Lal or I;(a*)=I; , for each j € J, we hereinafter focus on these 
columns. To do this, the following value matrix is defined based on problem 


(7)-(9). 


Definition 4. Define the value matrix of M=(mj;)mx2n, where 


c,(4;-«@;) if gi =1, 0 if g=1, 
™M4,2j-1 = i J 3) dij : and M425 — dij . 
io) otherwise, co otherwise, 


for each 1 € J and j € J. 


We will employ the branch-and-bound method with the jump-tracking 
technique to solve problem (1)—(2) using the value matrix M. Since #; and &,, 
for each 7 € J, cannot be selected simultaneously along a branch, we should 
modify the branch-and-bound method. We consider three modifications on 
this method similar to [1] as follows: 


1. If we choose #; (or #;) to branch from one node to another node, then 
we never use @; (or £;) to branch further on the current node. 
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2. Under the stated conditions below, we cannot branch further on Node 
k. 
2.1. We have reached to the last row of the matrix M. 
2.2. The selected variables along Node 0 to Node k together with £;, 
for each 7 € J \ Jz, satisfy all the equations, where 
Jy = {j € J\x; has been selected along the branches from Node 0 to Node k}. 
2.3. We do not have any candidate for satisfying an equation with 
regard to modification 1. 


3. If we cannot branch further on Node k under the conditions 2.1 and 
2.2, then we assign %; to x; for each j € J \ Jp. 


Note that if we cannot branch further on Node k with the value of Z, under 
the conditions 2.1 and 2.2, then Z, represents the objective value of problem 
(7)-(9) for the obtained vector x along Node 0 to Node k. Then the total of 
Zr along the branches from Node 0 to Node k can be calculated as follows: 


Total Z, = Z, + Ye; (10) 
jes 


where total Z;, represents the objective value of problem (1)—(2) for the ob- 
tained vector x. 

In problem (7)—(9), each equation i, for 7 € Jy, of its constraints can be 
satisfied by £;, for 7 € J; , or £;, for 7 € J;*. In this case, the ith equation 
may be satisfied without imposing any extra cost to the objective function. 
On the other hand, each equation i, for 7 € Jz, can only be satisfied by #;, 
for 7 € ie , that is, we have to expend extra cost for satisfying each equation 
i, for 7 € Ig. Since each equation i, for 7 € Ig, is satisfied with extra cost, 
we start the modified branch-and-bound method from row(s) 7 € Ig. In this 
case, the ith equation, for 7 € [,, may be satisfied with the expended cost for 
i € Ig or without imposing any extra cost to the objective function. Hence, 
we expect that the visited nodes of the modified branch-and-bound method 
can be decreased when we consider Remark 3. 

We are now ready to design an algorithm to solve problem (1)—(2) based 
on the obtained results up to now. 


4.2 An algorithm for the resolution of problem (1)—(2) 


Algorithm 1. Consider the optimization problem (1)-—(2). 

Step 1. Compute the lower and upper bound of # and # applying Lemma 
3. 

Step 2. If b=0 and #; < #;, for each j € J, then S(AT,A™,b) = {x |#< 
x < &} and x* = & is an optimal solution of problem (1)—(2) with regard to 
Corollary 1 and stop! 
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Step 3. If x; < %;, for each 7 € J, and 6; > 0, for each 7 € J, then go to 
Step 4. Otherwise, use Lemmas 4 and 6. 

Step 4. Compute the matrices of Q* and Q~, the index sets of iy and 
I; , for each j € J, and the index sets of J;* and J;, for each i € I using 
Definitions 1 and 2. 

Step 5. Check the consistency of bipolar max-T,, FREs (2) using Theorem 
2. If it is inconsistent, then stop! Otherwise, go to Step 6. 

Step 6. Perform the process of problem reduction as follows: 

6.1. Compute two index sets J; and J, using Definition 3. 

6.2. If If = 0, then «* = & is an optimal solution of problem (1)-(2) 
with regard to Corollary 3 and stop! 

6.3. If the conditions of Lemma 10 are satisfied, then the unique 
optimal solution x* of problem (1)—(2) can be obtained by relation (4) and 
stop! 

6.4. Check the conditions of Theorem 7. If the conditions are satisfied, 
then there exists an optimal solution x* according to relation (4) and stop! 

6.5. Check the conditions of Theorem 4. If the conditions are satisfied, 
then remove all the equations with numbers i € T;, from the matrices of QT 
and Q7 with regard to Corollary 2. 

6.6. If there exists a pair i € I and k € J such that g/,=q7,=1, then 
we can remove the row of 7 in the computation of the minimum objective 
value with regard to Theorem 6. 

Step 7. If Q‘=Q" =0, then assign %; to x} and go to Step 10. 

Step 8. Rearrange the rows of the matrices Qt and Q~ according to Remark 
3. Also, generate the value matrix of MW using Definition 4. 

Step 9. Employ the modified branch-and-bound method with the jump- 
tracking technique on the matrix M to solve the optimization problem of 
(1)-(2). 

Step 10. Produce the optimal solution and the optimal value of problem 
(1)-(2). End. 


5 Numerical examples 


We now illustrate Algorithm 1 by the following examples. 


Example 1. Consider the following optimization problem: 


min t1+ 3x2 + 273 + 5x4 + 845+ 726, (11) 
st. AtogxVA onr=b, (12) 
x € (0, 1", 


where 
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0.4 0.1 0.25 0.29 0.18 0.07 
0.32 0.16 0.23 0.1 0.2 0.48 
0.08 0.05 0.02 0.1 0.03 0.06 
0.15 0.3 0.17 0.2 0.14 0.05 | ’ 
0.09 0.2 0.03 0.04 0.15 0.24 
0.49 0.58 0.6 0.37 0.75 0.54 


Ae 


0.21 0.4 0.19 0.27 0.12 0.04 
0.15 0.32 0.8 0.14 0.22 0.2 
0.1 0.08 0.05 0.18 0.07 0.06 
0.17 0.04 0.15 0.07 0.3 0.2 
0.11 0.1 0.03 0.14 0.14 0.07 
0.24 0.52 0.47 0.5 0.33 0.25 


b = (0.3, 0.24, 0.09, 0.18, 0.12,0.6)7, and x = (%1,%2,23,24,%5,2%6)'. Now, 
we apply Algorithm 1 to solve the optimization problem of (11)-(12). 

Step 1. The lower and upper bound of & and & are as follows: 

# = (0:1,0.25, 0.7, 0:5, 0.4, 0.1)" and. % = (0:75, 06, 1,0,9,0.8, 0.5)". 

Step 2. Since the conditions of Corollary 1 do not hold, we go to Step 3. 
Step 3. In this example %; < #;, for each 7 € J and b; > 0, for each 7 € I. 
Therefore, we go to Step 4. 

Step 4. Applying Definition 1, the matrices of Qt and Q~ are obtained as 
follows: 


123456 123456 

1 /100000 1 /010000 
2{100001 2 {011000 
Qt = 3 | 000100} andQ- = 3] 100100 
41010100 41000011 

5 |010011 5 | 000000 

6 \001010 6 \000000 


Also, the index sets of iy and I; , for all 7 € J can be computed as 
follows: 
If = {1,2}, It = {4,5}, Ix = {6}, It = {3,4}, It = {5,6}, Ig = 
{2,5}, Ty = {3}, 17 = {1,2}, Iz = {2}, Ig = {3}, i {4}, and 
Ig = {4}. 

Moreover, we can compute the index sets of J;* and J;, for all i € J, as 
follows: 
Jt = {1}, Jy = {1,6}, Jz _ 44}, ay = {2,4}, Js _ {2, 5, 6}, Ig = 
{o,f S12}, do S125) d, ={1h4), 47 = {6,6),,4nd JC ade =o. 
Step 5. The bipolar max-T, FREs of At oz V A~ o72@ = b is consistent 
according to Theorem 2. So, we go to Step 6. 
Step 6. Perform the process of problem reduction as follows: 
6.1. Two index sets Jy and Ip are as follows: I; = {1,2,3,4} and In = {5,6}. 
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Since the optimization problem of (11)—(12) cannot be reduced by Steps 
6.2-6.4, we go to Step 6.5. 
6.5. In this example, T; = {2}, that is, Jj C Jf and J; C Jz. Applying 
Corollary 2, the second row of two matrices Q* and Q~ can be removed. 
6.6. Since aa, = q34 = 1, the third equation can be eliminated from our 
consideration with regard to Theorem 6. 

The matrices of Q* and Q~ cannot be reduced further. So, we go to Step 
ve 
Step 7. Since Qt 4 and Q~ £9, we go to Step 8. 
Step 8. With regard to Remark 3, the matrices of Qt and Q~ can be 
updated as follows: 


123456 123456 
5 (010011 5 (000000 
Qt+=6]001010)andQ- = 6 | 000000 
1 | 100000 1 | 010000 
4 \010100 4 \000011 


For the updated matrices of QT and Q~, the value matrix of M can be 
generated as follows: 


1 2 3 4 5 6 
co co 61.050 wWMmHo wo 3.20 2.80 
M= CO CO co CO (0600 WO O82 


0.65 co co 0 COCO WKH WW Wo 
co Cc) «61.0500 wo 2c ow 0 co 0 


Step 9. We are now ready to use the modified branch-and-bound method 
with the jump-tracking technique on the matrix M. We begin with the first 
equation, that is, i = 1. The set of {@2,%5,%¢} introduces three candidates 
to satisfy the first equation. Therefore, we have to branch from Node 0 in 
Figure 2. If we select #2: (Node 1), then the value of Z is 1.05. Note that 
%2 cannot be used for further branching on Node 1. Also, if we select #5 
(%6), then we have Z2 = 3.2 (Z3 = 2.8). Furthermore, we never use #5 (Ze) 
to branch further on Node 2 (Node 3). Here, Node 1 is selected to branch 
further because of the least objective value. 

Move to the second row of the matrix M. Since the set of {%3, #5} contains 
two candidates to satisfy the second equation, we have two branches from 
Node 1 as it has been illustrated in Figure 2. If #3 (Node 4) is selected, then 
the value of Z4 is 1.65. If #5; (Node 5) is considered, then we have Z5 = 4.25. 
Now, we can branch further on four Nodes 2, 3, 4, and 5 but Node 4 is chosen 
with regard to the least objective value. 

Move to the third row of the matrix M. The set of {%1, %2} contains two 
candidates to satisfy the third equation, but %2 cannot be used to branch 
further on Node 4 because <2 has been chosen along Node 0 to Node 4 
(modification 1). Therefore, the set of {#1} contains the only candidate to 
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Figure 2: The modified branch-and-bound method 


satisfy the third equation here. If #, (Node 6) is selected, then it is concluded 
that Z = 2.3. 


Since £1, £2, and <3 together with #4, @;, and Z&@ satisfy all the equa- 
tions, we do not branch further on Node 6 with regard to modification 2.2. 
Therefore, considering modification 3, the vector r=(#1, £2, £3, £4, £5, %6)" 
is a solution with the objective value of 2.3 for the equivalent problem. Since 
the value of Ze is less than Z2, Z3, and Zs, we can stop further branching 
on all Nodes 2, 3, and 5. Thus, the total of Z can be computed as follows: 

6 


Total Z6=Ze + Y> cj%; = 2.3 + 8.65 = 10.95. Every node is stopped further 
j=l 


branching. Hence, the optimal solution can be obtained from Node 6, that 
is, xj = 0.75, 5 = 0.6, 3 = 1, 7] = 0.5, e§ = 0.4, and x§ = 0.1 with the 
total of Z = 10.95. 

Step 10. The optimal objective value is 10.95 and the optimal solution is 
eo = (ot eel eek) = (075, 06,1,05,04,01)". 

If we do not apply Remark 3 to solve this example, we need to consider 25 
nodes. Considering Remark 3, we need only six nodes to solve this example. 


Example 2. Consider the following optimization problem: 


min 22, + 5x2 + 3%3 + 404+ 25 + 6246, (13) 
st. AtogxVA onr=b, (14) 
we (0, 1)", 


where 
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0.7 0.3 0.96 0.66 0.8 0.05 
0.31 0.54 0.34 0.3 0.36 0.45 
0.06 1 0.04 0.625 0.03 0.08 
0.15 0.15 0.19 0.225 0.16 0.03 | ’ 
0.6 0.08 0.07 0.06 0.13 0.4 
0.3 0.05 0.09 0.07 0.11 0.19 


Ars 


0.11 0.4 0.19 0.27 0.12 0.04 
0.14 0.24 0.09 0.14 0.22 0.2 
0.1 0.08 0.05 0.18 0.8 0.8 
0.2 0.04 0.15 0.3 0.11 0.05 | ’ 
0.11 0.4 0.32 0.14 0.14 0.07 
0.03 0.2 0.06 0.02 0.04 0.07 


b = (0.6, 0.27, 0.5, 0.18, 0.24, 0.12)7, and x = (x1, 22,23, 24,%5,2%6)'. Now, 
we apply Algorithm 1 to solve the optimization problem of (13)-(14). 

Step 1. The lower and upper bound of & and & are as follows: 

# = (0:1, 0.4, 0.25, 0.4, 0.375, 0.375)" and #= (0.4, 0.5, 0.625, 0.8, 0.75, 0.6)7. 
Step 2. Since the conditions of Corollary 1 do not hold, we go to Step 3. 
Step 3. In this example, £; < £;, for each 7 € J and b; > 0, for each 7 € I. 
Therefore, we go to Step 4. 

Step 4. Applying Definition 1, the matrices of Qt and Q~ are obtained as 
follows: 


123456 123456 

1 /001010 1 /000000 
21010011 2 {000000 
Qt = 3 |010100]andQ- = 3] 000011 
41000100 41100100 

5 |100001 5 |011000 

6 \100000 6 \010000 


Also, the index sets of if and I; , for all 7 € J can be computed as 
follows: 
TY = {5,6}, 1 = {2,3}, B = (}, Wt = {3,4}, TP = {1,2}, = 
{2,5}, = {4}, Ig = {5,6}, Ig = {5}, Ib = {4}, Is = {3}, and 
Ig = {3}. 

Moreover, we can compute the index sets of J;* and J;, for all i € J, as 
follows: 
FF = {3,5}, JY = {2,5,6}, Js = {2,4}, Jd = {4}, JP = {1,6}, Je = 
{1}, J, =Jz =0, Jz = {5,6}, Jg = {1,4}, Js = {2,3}, and Jg = {2}. 
Step 5. The bipolar max-T, FREs of At oz V A~ o7@ = b is consistent 
according to Theorem 2. So, we go to Step 6. 
Step 6. Perform the process of problem reduction as follows: 
6.1. Two index sets I, and Jy are as follows: I, = {3,4,5,6} and Ip = {1,2}. 
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6.2. The condition of Substep 2 is not satisfied for this problem. 

6.3. The conditions of this substep hold for this problem. Since for each 
9 6d =e Slt Of — (Ji U Jt) = {1,4}, we have c; = 2 > 0 and 
c4 =4 > 0. Also, for k € (),<7, = {5}, the following conditions are satisfied: 
1. forall j € Use; Jf — {5} = {2,3,6}, we have c5(#5 — #5) = 0.225 < 
cj(#j — #5). 

2. Is — Is = {3} CUjerj¢5 1; = {3,4,5; 6}. 

Therefore, according to Lemma 10, the unique optimal solution of the prob- 
lem is as follows: z* = (0.1, 0.4, 0.25, 0.4, 0.75, 0.375) with the optimal objec- 
tive value Z* = 7.55 


6 Comparison of the proposed algorithm with the 
methods in other papers 


In this section, we compare the proposed algorithm with the methods in other 
papers [9, 20, 23, 3] to solve problem (1)—(2) with regard to the obtained 
results from Examples 1 and 2 in Section 5. 

As it is seen in Figure 2, Algorithm 1 solves the problem of Example 1 
only in six nodes by the branch-and-bound method. Also, Algorithm 1 solves 
Example 2 without using the branch-and-bound method. Its optimal solution 
is found in Substep 6.3 by Lemma 10. 

Freson, De Baets, and De Meyer [9] discussed problem (1)—(2) with the 
max-min composition operator. They designed an algorithm to solve the 
problem. Its Step 1 checks the necessary condition (46) in [9]. Its compu- 
tational cost is O(mn). Step 2 constructs vectors gt and s~ taking supre- 
mum and infimum from m maximal solutions which its computational cost 
is O(mn). Step 3 generates all elements of set B and keeps those that satisfy 
constraints (30) in [9]. To generate all the elements of B according to relation 
(48) in [9], we firstly need to produce the set ({0, 1} U{b;, b;|4 = 1,...,m})”. 
Its computational cost is O((2m+2)”). Then we must check whether each its 
element belongs to [s~,g*] or not? Its computational cost is O(n(2m +2)"). 
Now, we must check whether each element B satisfies constraints (30) in 
[9] or not? The computational cost of this work is O(mmn) for each ele- 
ment of the set B. If we check all the elements of set B, then its com- 
putational cost is O(mn(2m + 2)"). So, the computational cost of Step 3 is 
T3 = O(n(2m+2)"+mn(2m+2)”) = O(mn(2m+2)”). Step 4 selects the el- 
ements in B(| D with the highest value for the objective function. To do this, 
we must check the objective function for |B (| D| elements which its computa- 
tional cost is O(|B() D|) < O((2m+2)”). Therefore, the computational com- 
plexity of the given algorithm in [9] is Tr = O(mn(2m+2)”). For an instance 
of the problem with the dimensions m = n = 6 like Examples 1 and 2, its 
computational complexity is Tr = O(6 x 6 x (2x 6+2)°) = 271063296 x O(1). 
This point implies that (2 x 6 + 2)© = 7529536 elements are produced and 
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its feasibility is checked in m equations of n-variable. Then the optimizer is 
found by computing the objective function values in the feasible vectors. 


In [20], problem (1)—(2) is discussed with the max-Lukasiewicz t-norm 
composition. Li and Liu [20] directly converted the problem to a 0-1 integer 
linear optimization problem without its simplification and reduction. If we 
apply their method to solve the problem of Example 1, we should consider 
the following 0-1 integer programming problem: 


Z=min 8.65+0.65u; + 1.05ug + 0.6u3 + 2u4 + 3.2u5 + 2.8u6, 


tO 0 0\ Fe 0 
i <i oa 4) |e =i 
= 0 GO: 6 || | as = 

as 01 01-1-1] |u| =]—1]? 
oe ae eee ee 1 
Oo 10d Oe 1 
u € {0,1} 


If we apply Algorithm 1 of this paper for Example 1, we should solve the 
following problem with smaller dimensions: 


Z=min 8.65+0.65u; + 1.05ug + 0.6u3 + 2u4 + 3.2u5 + 2.8u6, 


OL00 tT T\ [uw 1 
00101 0 | Jus 1 
> 
St 1119000 o}]lu|=lo |? 
1.01 a, =i 


u € {0,1}°. 


If we apply Li and Liu’s method to solve the problem of Example 2, we should 
consider the following 0-1 integer programming problem: 


Z=min 7.175 + 0.6uz + 0.5u2 + 1.125uzg + 1.6u4 + 0.375u5 + 1.35u6, 
004101 #0 U1 1 
01001 1 ug 1 
0 1 01-1-1 U3 —1 
oe lait 00 @@ | au | i$ 
1-1-1000 1 Us —1 
11-1000 0 U6 0 
u € {0,1}°. 


If we apply Algorithm 1 of this paper for Example 2, we directly obtain 
the unique optimal solution of the problem of Example 2 without using the 
branch-and-bound method and 0-1 integer programming problem. 
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In [3], problem (1)—(2) was discussed with the max-Hamacher t-norm com- 
position with some rules to simplify the problem. The rules are completely 
different from the proposed procedures of this paper for simplification. The 
rules in [3] are not applicable for the problem in Example 1. Hence, if we 
use the given algorithm in [3] for Example 1, the branch-and-bound method 
should be applied on the matrix M as follows: 


1 2 3 4 5 6 
0.65 co co 0 O00, CEEC O10 © O10 0REEEC 0110.6) 
0.65 co coo 0 o0 ww wo. 2.8 co 
M= oo 0 0 0 6K Oo 0 OC (15) 
co co 61.0500 wo 20 ow 0 coo 0 | 
co co «61.0500 Wom wo 3.20 2.80 
CO OO co co 60600 WC 8.20 


If we use Algorithm 1 for Example 1, the branch-and-bound method should 
be applied on the matrix MW as follows: 


1 2 3 4 5 6 
co cf0 §61.0500 wo Woo 3.20 2.800 
M= 00 00 —60 CON.B COO Hs 38.2 WW” (16) 


0.65 co co 0 (cone OREO 00 OREO Ole ©MEEE® 0110 © 
co co) «61.0050 wo 20 ow 0 co 0 


Also, the rules in [3] are not applicable for the problem in Example 2. Hence, 
if we use the given algorithm in [3] for Example 2, the branch-and-bound 
method should be applied on the matrix MW as follows: 


1 2 3 4 5 6 
coo 6m oo) 1.12500 wo 0.3875 co oO CO 
co co) 60.5 co co co wc. 0.3875 0c 1.35 co 
M= | wo 0.500 co co 1.600 oo «(0 oo 0 (17) 
co 0 co CO co coo 1.60 CO CO co Co | 
0.60 ow 0 co 0 oO OO co co 1.3500 
0.60 ow 0 co Cc «co co CO oO CO 


If we use Algorithm 1 for Example 2, its optimal solution is directly found in 
Substep 6.3 without using the branch-and-bound method and M = 9. 


In [23], problem (1)—(2) was discussed with the max-Lukasiewicz com- 
position with some rules to simplify the problem. The rules are completely 
different from the proposed procedures of this paper for simplification. If we 
apply the method in [23] for the problem of Example 1 with the max-product 
composition, rule 3 in [23] can be used for its simplification. Applying the 
equivalent form of the rule 3 for problem (1)-(2), the forth row of the matrix 
M in the relation (15) is removed and the branch-and-bound method should 
be employed on the 0-1 integer programming equivalent to the following value 
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matrix: 
1 2 3 4 5 6 
0.65 oo oo oe EEC OR OMEN O10 0G O10 ©) 
we 0.65 co oo co 0 wo wo.) 2.8 0c (18) 
co (0 oo 


C0 CO 6B 8.200 2.8 CO 


0 
0 
o wc 2 0 ooo oO Oe ©) 
oe) 
o 060 wo 3.20 Wo 


If we use Algorithm 1 for Example 1, the branch-and-bound method should 
be applied on the matrix M in the relation (16) and considering Remark 3, 
we need only six nodes to solve this example. If we apply the method in 
[23] for the problem of Example 2 with the max-product composition, rule 
3 in [23] can be used for its simplification. Applying the equivalent form 
of rule 3 in [23], the third row of the matrix M in relation (17) is removed 
and the branch-and-bound method should be employed on the 0-1 integer 
programming equivalent to the following value matrix: 


1 2 3 4 5 6 
coo «6M 1.12500 = ow O0..875 co oO CO 
aps co cCO-)=6—(0.5 co oOo co Mc 0.3750 1.3500 (19) 
co 0 oO OO co ow 1.6 0 oO CO oO CO 
0.6c ow 0 co (0 oO OO co ow 1.35 00 
0.600 ow 0 co CO Wc CO 60CO co CO 


If we use Algorithm 1 for Example 2, its optimal solution is directly found in 
Substep 6.3 without using the branch-and-bound method and M = 9. 


The other preferences of the proposed algorithm with respect to the pre- 
sented algorithms in [9, 20, 23, 3] for the resolution of problem (1)—(2) are 
as follows. The proposed algorithm introduces two new classes of problem 
(1)-(2) with the max-product composition operator, which can directly be 
solved only by Substep 6.3 or Substep 6.4. The optimal solution of these 
classes of the problem satisfying conditions of Substep 6.3 or Substep 6.4 
can be obtained by the relation (4) without applying the branch-and-bound 
method or using rules repeatedly. The classes have not been introduced in 
(9, 20, 23, 3]. Other proposed rules in Step 6 of Algorithm 1 are different 
from the given rules in [23, 3]. In [9, 20], the authors have not used the 
rules of simplification to reduce the original problem. Since the algorithms 
in [20, 23, 3] and the proposed algorithm are based on the branch-and-bound 
algorithm, the algorithms are convergent. The algorithm in [9] checks all the 
possible feasible solutions to find the optimal solution. 
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7 Conclusions and future works 


The linear optimization problem with the bipolar max-product FREs was 
studied in this paper. The characterizations of its feasible domain were in- 
vestigated. Some simplification operations were proposed to delete some 
equations. With regard to these operations, the size of the original problem 
was reduced. Then, some sufficient conditions were presented to determine 
one of the optimal solutions to the problem and its uniqueness. Moreover, a 
value matrix was defined based on the characteristic matrices of the feasible 
domain of the problem. Then, the branch-and-bound method was modified 
to solve the reduced problem with regard to the value matrix. An algorithm 
was finally designed to solve the problem and compared with other methods 
to show the efficiency of the proposed algorithm. In future work, the linear 
optimization problem will be developed by supposing fuzzy linear systems 
with bipolar fuzzy numbers based on references [10, 11]. 
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